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The charging of a quantum box connected to a lead by a single-mode point contact is solved 
for arbitrary temperatures, tunneling amplitudes, and gate voltages, using a variant of Wilson's 
numerical renormalization group. The charge inside the box and the capacitance of the junction 
are calculated on equal footing for all physical regimes, including weak tunneling, near perfect 
transmission, and the crossover regime in between. At the charge plateaus, perturbation theory is 
found to break down at fairly small tunneling amplitudes. Near perfect transmission, we confirm 
Matveev's scenario for the smearing of the Coulomb-blockade staircase. A surprising reentrance of 
the Coulomb-blockade staircase is found for large tunneling amplitudes. At the degeneracy points, 
we obtain two-channel Kondo behavior directly from the Coulomb-blockade Hamiltonian, without 
the restriction to two charge configurations or the introduction of an effective cutoff. 

PACS numbers: 73.23.Hk, 72.15.Qm, 73.40.Gk 



The Coulomb blockade^ is one of the fundamental 
phenomena in mesoscopic physics. When a quantum box, 
either a small metallic grain or a large semiconducting 
quantum dot, is coupled by weak tunneling to a lead, its 
charging is governed by the finite energy barrier Ec for 
adding a single electron to the box. This gives rise to the 
well-known Coulomb-blockade staircase for the charge of 
the box as a function of gate voltage^ Increasing the 
coupling to the lead smears the Coulomb staircase, as 
thermal fluctuations do at ksT > Ec- 

A large diversity of theoretical approaches have been 
applied to the Coulomb blockade, ranging from pertur- 
bation theory&i^ and diagrammatic techniques to 
renormalization-group treatments^iiSiii*^ and Monte 
Carlo simulations>ii*i2ii2ii^ However, to date there is no 
single unified approach encompassing all regimes of the 
Coulomb blockade. The reason being that different start- 
ing points are required for describing the different physi- 
cal regimes of the Coulomb blockade. For weak tunneling 
and UbT <C Ec, the charge inside the box is essentially 
quantized at the charge plateaus. Hence, one can start 
from a well-defined charge configuration and apply the 
perturbation theory in the tunneling amplitude. This 
approach collapses at the degeneracy points between ad- 
jacent charge plateaus, where strong charge fluctuations 
give rise to exotic many-body physics in the form of the 
two-channel Kondo effecti£*i£ For large transmission, the 
Coulomb-blockade staircase is smeared out and the no- 
tion of well-defined charge configurations breaks down. 

In this paper, we devise such a unified approach for 
all temperatures and parameter regimes of the Coulomb 
blockade using Wilson's numerical renormalization-group 
method. 1 - Focusing on single- mode point contacts, we ac- 
curately calculate the charge and the capacitance of the 
quantum box, for arbitrary gate voltages, temperatures, 
and tunneling amplitudes. These quantities were recently 
measured by Berman et al. for single-mode point con- 
tacts^ revealing signatures of the two-channel Kondo 
effect. For weak tunneling, we recover the results of the 
perturbation theory at the charge plateaus. However, we 



find that the perturbation theory breaks down at fairly 
small tunneling amplitudes. Close to perfect transmis- 
sion, we find good agreement with Matveev's analysis of 
a related model^ which in turn breaks down as one de- 
parts from perfect transmission. Upon further increasing 
the tunneling matrix element, there is a surprising reen- 
trance of the Coulomb-blockade staircase. At the degen- 
eracy points, we obtain two-channel Kondo behavior di- 
rectly from the Coulomb-blockade Hamiltonian, without 
the restriction to two charge configurations or the intro- 
duction of an effective cutoff. Indeed, the effective cut- 
off entering Matveev's two-charge-configuration model? 
is found to depend not only on the charging energy Ec, 
but also on the tunneling strength. 

The coupled quantum-box-lead system is convention- 
ally modeled by the Hamiltonian 

H = E E C «*4*aCa*a+ ^ + V B Qb (1) 
a — L,B ka 

+ 1 E { C lka C Bk>* + H.C.) , 
kk' a 

where c* Lka (c Bka ) creates a lead (box) electron with mo- 
mentum k and spin a; t is the tunneling matrix element, 
taken for simplicity to be momentum independent; and 
tLk (esfc) are the single-particle levels in the lead (box), 
which are assumed to be continuous (dense energy lev- 
els). The excess electrical charge inside the box, 

Qb = -e Yy*Bka c Bka - 0{-e Bk )] (2) 

ka 

(— e being the electron charge), is controlled by the ca- 
pacitance of the box, Co, and by the external component 
of the electrostatic potential in the box, V B ■ 

The numerical renormalization groupie (NRG) is a 
nonperturbative approach for treating quantum impurity 
problems. At the center of this approach is a logarith- 
mic energy discretization of the conduction band around 
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N B = eV B /2E c 



FIG. 1: The excess charge inside the box versus Vb, for 
tp — 0.1 and different temperatures. Here, Ec/D = 0.01, 
A = 2.5, and N s = 2000. At temperatures comparable to 
the charging energy, the charge curve is essentially linear. A 
sharp Coulomb-blockade staircase emerges for fcsT <C Ec- 
Inset: The capacitance versus Vb, for the same temperatures 
as plotted in the main panel. For ksT > Ec, the capacitance 
is constant and equal to Co = e 2 /2Ec- For fcsT <S Ec, sharp 
capacitance peaks develop at the degeneracy points. 



the Fermi level, designed to capture all logarithmic diver- 
gences of the problem. Using a unitary transformation, 
the conduction band is mapped onto a semi-infinite chain 
with the impurity coupled to the open end, and a hop- 
ping matrix element that decreases exponentially along 
the chain. A direct application of this approach to the 
Coulomb-blockade Hamiltonian is doomed to fail, since 
the Coulomb interaction equally couples all sites along 
the semi-infinite chain. It does not decay as the hopping 
matrix elements. To circumvent this problem, we resort 
to a mapping used by Schoeller and co-workers Ex- 
plicitly, we map the Hamiltonian of Eq. Q onto 



H = 



e akC ! akcr C a kcr 

x—L.B ka 



t 



E 

kk' a 



j 

^Lka 



CBk>, 



2C 
H.c.}, 



V B Q (3) 



where 



\n±l)(n\ and Q — — e^^n\n){n\ (4) 



are new collective charge operators. Strictly speaking, 
the above mapping requires the constraint Q = Q B . 
However, this constraint can be relaxed in the contin- 
uum limit, when the dynamics of Q in Eq. @ is insensi- 
tive to the precise number of conduction electrons in the 
bands. Hereafter, we regard Q as an independent degree 
of freedom. 

Equation describes two noninteracting conduction 
bands, coupled to an impurity with an infinite number 



of parabolically dispersed energy levels. Hence, it can 
be treated using the NRG. We solved this model using a 
common constant density of states p for the lead and the 
boxjiS and a common bandwidth D. The NRG method 
is controlled by the logarithmic discretization parameter 
A and by the number of states retained N s m& 

Figure n shows the expected Coulomb-blockade stair- 
case at different temperatures, for the weak tunneling 
tp = 0.1. At temperatures above the charging energy 
Ec = e 2 /(2Co), one recovers a linear (classical) charge 
curve. As the temperature is reduced below the charg- 
ing energy, a staircase structure emerges, with a period 
equal to twice the charging energy. For fc^T <C Ec, 
one is left with sharp charge steps, with abrupt transi- 
tions between the charge plateaus at half-integer values 
of N B = eV B /2E C - 

A complimentary picture is provided by the capaci- 
tance, C = —d{Q)/dVB, recently measured by Berman 
et alJL As seen in the inset of Fig. ^ the capacitance C 
is equal to the classical capacitance Co, for k B T > Ec- 
However, C acquires a strong dependence on N B , for 
k B T <C Ec- While C is significantly reduced at the 
charge plateaus, near the degeneracy points (i.e., near 
half-integer values of N B ) it is dramatically enhanced. 
As will be discussed below, C diverges logarithmically at 
the degeneracy points as T — * 0, in accordance with the 
two-channel Kondo effect predicted by Matveevi^ 

At the charge plateaus, perturbation theory in tp is 
commonly used for weak tunneling. To second order in 
tp, the zero-temperature excess number of electrons in- 
side the box, n = — (Q)/e, is given by 



n = 2(tp) 2 \ In 



1 + 2N B 



2 In 



1-2A^b 
1 - 2N B + d 
1 + 2Af B + d 



(5) 



In 



1 + 2N B + 2d 
1 - 2N B + 2d 



Here, we have considered particle-hole symmetric bands, 
and retained a finite ratio d = D / Ec- In the wide-band 
limit, d — > oo, the last two terms drop in Eq. JSJ, and one 
recovers the conventional expression (see, e.g., Ref. 0). 

The inset of Fig.[3compares between the NRG and the 
perturbative expression of Eq. JSJ . At very weak tunnel- 
ing, tp = 0.05 and 0.1, the NRG results coincide with the 
perturbative expression. The good agreement extends to 
all gate voltages except the vicinity of the degeneracy 
points. However, significant deviations from Eq. J^J are 
seen at fairly weak coupling, tp w 0.2 — 0.3. Indeed, the 
Coulomb-blockade staircase is completely washed out for 
tp = 0.3 (main panel of Fig.[2J), indicating the breakdown 
of the perturbation theory. Moreover, there is a reen- 
trance of the Coulomb-blockade staircase for tp > 0.4, 
which obviously is not captured by the perturbation the- 
ory in tp. 

The rapid breakdown of the perturbation theory is 
related to the approach to perfect transmission. For 
Ec = and particle-hole symmetric bands, the single- 
particle transmission coefficient at the Fermi level is given 
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FIG. 2: Evolution of the zero-temperature charging curve 
with increasing tp. Here, Ec/D — 0.01, A = 2.5, and 
N s — 2000. The values of tp are specified in the legend. 
At very weak coupling, there is a sharp Coulomb-blockade 
staircase. As demonstrated in the inset, good agreement is 
obtained with the perturbative expression for n/2(tp) 2 [see 
Eq. J5J], which is plotted for comparison by the solid line. 
Significant deviations from the perturbation theory develop 
for tp > 0.2. Specifically, the Coulomb-blockade staircase is 
completely washed out for tp — 0.3, marking the breakdown 
of the perturbation theory. For tp > 0.4, there is a reentrance 
of the Coulomb-blockade staircase. 



by 



T = 4- 



{ntpf 



[1 + (irtpy 



(6) 



For ntp = 1, T reaches perfect transmission and is no 
longer perturbative in tp. That T rather than tp is the 
relevant physical parameter for Ec <C D is evident from 
Fig. The Coulomb- blockade staircase is washed out for 
tp ~ 0.3 — 0.4, when T of Eq. (JSJ approaches unity. The 
reentrance of the Coulomb-blockade staircase for tp > 0.4 
stems from the reduction of the single-particle transmis- 
sion coefficient with increasing tp > 1/vr. It should be 
emphasized, however, that this picture is only valid for 
E c <C D. Upon increasing E c , the charge curves can- 
not be exclusively parametrized by the transmission co- 
efficient of Eq. (JBJ. In particular, the washing out of 
the Coulomb-blockade staircase and its subsequent reen- 
trance are pushed to higher values of tp, reflecting the 
limitations of the single-particle picture. 

Near perfect transmission, it was predicted by Matveev 
that the zero-temperature excess number of electrons in- 
side the box is given by 1 - 



n = N 



b + - J1 j~ In [e\r\ 2 cos 2 {ttN b )] sin(2^V B ) . (7) 



Here, \r\ 2 = 1 — T is the reflectance, 7 is equal to e c , and 
C k, 0.5772 is Euler's constant. This result was obtained 
within a one-dimensional model, especially designed for 
T w 1 . This regime of large transmission was never stud- 
ied using a tunneling Hamiltonian. 




FIG. 3: Comparison of the NRG (symbols) and Eq. Q, for 
the zero-temperature charge curve near perfect transmission. 
Here, E c / D = 0.01, A = 2.5, and N B = 2000. The horizontal 
axis, Sn — n — Nb, measures the deviation from linearity 
of the excess number of electrons inside the box. The solid 
lines are fits to Eq. Q, using the single fitting parameter \r\ 2 . 
Inset: A plot of the extracted values of \r\ 2 versus tp. The 
solid line shows 1 — T of Eq. J5J with tp — > tp — 0.023. 



Figure [3] shows a comparison of our NRG results for 
Sn = n — Nb and Eq. J7J), using the single fitting pa- 
rameter |r| 2 . There is a good agreement between the 
two curves, particularly for 0.3 < tp < 0.4. Moreover, 
the extracted values of \r\ 2 closely trace 1 — T in this 
range of tp (see inset of Fig. |3J), apart from a slight shift 
which is likely due to the finite value of Ec / D used. For 
tp = 0.25 and tp — 0.45, small systematic deviations 
from Eq. (JJJ develop in the middle of the plateaus, near 
integer values of Nb- Further away from perfect trans- 
mission, e.g., for tp = 0.2 and tp — 0.5 (not shown), there 
are large deviations from the analytic expression, signal- 
ing the breakdown of Matveev's expression. This clearly 
shows that one can still use the tunneling Hamiltonian 
of Eq. near perfect transmission, provided that the 
tunneling processes are summed to all orders. 

Of particular interest is the behavior near the degener- 
acy points, where the perturbation theory diverges log- 
arithmically even for small values of tp. For weak tun- 
neling and ksT <C Ec, the charge fluctuations in this 
regime were mapped by Matvee\£ onto the planar two- 
channel Kondo model, by retaining only the two lowest- 
lying charge configurations inside the box. The capaci- 
tance at the degeneracy points was predicted to diverge 
logarithmically with T — > 0, in accordance with the log- 
arithmic divergence of the magnetic susceptibility in the 
two-channel Kondo effect. 15 This scenario was never con- 
firmed directly for the Hamiltonian of Eq. (JIJ. 

Figure 0] shows the capacitance at the degeneracy 
points as a function of temperature. At low tempera- 
tures, the capacitance diverges logarithmically with de- 
creasing T, confirming the onset of the two-channel 
Kondo effect. Indeed, the two-channel Kondo effect per- 
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FIG. 4: The capacitance at the degeneracy points as a func- 
tion of temperature, for Ec/D — 0.25 and different values 
of tp. Here, A = 2 and N s = 2000. All curves show the 
characteristic logarithmic temperature dependence of the two- 
channel Kondo effect. Using Eq. @, we extract the Kondo 
temperature Tk from the slope of the logarithmically diverg- 
ing component of C(T). Inset: Tk/T k versus (tp) -1 , for 
Ec/D = 0.01 (open squares), 0.25 (open triangles), and 1 
(dots). Here, T K is the Kondo temperature for tp — 0.1. Ex- 
plicitly, k B T K /D equals 3.2 xlO" 6 , 3.6xl0~ 5 , and6.5xl0~ 5 , 
for Ec/D — 0.01, 0.25, and 1, respectively. 



where D c g is an effective high-energy cutoff, assumed to 
be of the order of the charging energy. 



In the inset of Fig.Q] we plotted T K /T K versus (tp) -1 , 
for a wide range of charging energies Ec/D. Here, T K 
is the Kondo temperature for tp = 0.1, extracted sep- 
arately for each value of Ec/D. For weak tunneling, 
Tk decays exponentially with (tp) , in accordance with 
Eq. JjJJ. In particular, all curves fall on top of one an- 
other in this regime, confirming that Ec enters the ex- 
pression for Tk only through the preexponential factor 
D c g. However, a detailed comparison of the NRG data 
and Eq. © reveals an additional dependence of _D e ff on 
tp. For Ec/D = 0.01, for example, D e ff is enhanced by 
a factor of 4 upon going from tp = 0.05 to tp = 0.2. 
At perfect transmission, the two-channel Kondo effect 
breaks down, which is signaled by a cusp in Tk as a 
function of tp [the slope of the ln(T) component of C(T) 
vanishes]. In contrast to the universality of the weak- 
tunneling regime, the position of the cusp in Tk shifts to 
higher values of tp with increasing Ec , showing that one 
cannot exclusively parametrize the system by the nonin- 
teracting transmission coefficient of Eq. JfjJ). 



sists for all values of Ec/D and all values of tp up to 
perfect transmission, where it breaks down^ However, 
the logarithmic behavior is regained upon further increas- 
ing tp, marking the reentrance of the two-channel Kondo 
effect. 

The crossover from the high-temperature to the low- 
temperature universal regime is generally marked by the 
Kondo temperature Tk, which we extract from the Bcthc 
ansatz expression for the slope of the logarithmically di- 
verging term in the capacitance: 2 ^ 



C(T) 



20k B T K 



HTk/T) . 



(8) 



Within the two-charge-configuration model of Matveev, 
Tk is given by^ 



T K = D cS tp exp 



Atp 



(9) 



In summary, we devised an NRG approach for solving 
the charging of a quantum box connected to a lead by 
a single-mode point contact, which uniformly treats all 
regimes of the Coulomb blockade. Using this approach 
we are able to (i) reveal the rapid breakdown of the per- 
turbation theory, (ii) obtain a surprising reentrance of 
the Coulomb-blockade staircase for large tunneling am- 
plitudes, (iii) confirm Matveev's scenario for the shape 
of the charge curve near perfect transmission, and (iv) 
obtain two-channel Kondo behavior at the degeneracy 
points directly from the Hamiltonian of Eq. • 
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